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Abstract 

We have calculated the spatial distribution of the electrostatic potential cre- 
ated by an unbalanced charge q in one of the conducting electrodes of a long, 
uniform, linear array of small tunnel junctions. The distribution describes, 
in particular, the shape of a topological single-electron soliton in such an 
array. An analytical solution obtained for a circular cross section model is 
compared with results of geometrical modeling of a more realistic structure 
with square cross section. These solutions are very close to one another, and 
can be reasonably approximated by a simple phenomenological expression. 
In contrast to the previously accepted exponential approximation, the new 
result describes the crossover between the linear change of the potential near 
the center of the soliton to the unscreened Coulomb potential far from the 
center, with an unexpected "hump" near the crossover point. 
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Recent theoretical and experimental studies have resulted in considerable progress in 
understanding correlated single-electron transfer in ultra-small tunnel junctions (for reviews, 
see Refs. These phenomena may be used as a background for a new generation of 

analog and digital devices. The most common component of single-electronic devices is 
a one- dimensional array of small tunnel junctions (Fig. [Ij(a) ) . Thus it is very important 
to gain a quantitative understanding of the Coulomb interaction potential, U(r), between 
single electrons in such an array. 

To our knowledge, all previous works on this topic (see, e.g., Ref. § and references therein) 
have used a simple model in which the complete matrix [C] of mutual capacitances between 
conducting "islands" of the array is truncated to tridiagonal form. In this form of the matrix, 
the only non-zero elements are (a) the diagonal elements = C , representing the stray 
capacitances of the islands, and (b) the nearest-neighbor elements C^i±i = C, dominated 
by tunnel junction capacitances. Electron-electron interaction in the tridiagonal model is 
described by a simple exponential lawSi: 

U t (r) = U t (0)exp(-m/m o ), (1) 

where m = r/a is the distance between the two electrons, in units of the array period a (i.e., 
in number of islands). The parameters U t (0) and m depend on the C/C Q ratio, and in the 
most important limit of C Q <C C: 

m ° = vf- (3) 

The tridiagonal model is strictly correct only if the array is placed parallel to, and very 
near, a conducting ground plane. However, we are not aware of any experiments which 
actually have used such a configuration. The presence of a ground plane would increase the 
stray capacitance C Q , thus cutting off the single-electron soliton radius m and suppressing 
the electron-electron interaction at large distances. Thus, it was our goal to describe long 
linear arrays of realistic geometry without a ground plane. 



If the length scale of the electron-electron interaction within the array is much larger 
than the array period (i.e., if m ^> 1), the electrostatics of the array should not depend 
strongly on the details of the geometry of its islands. Thus we can model the experimental 
array geometry, such as the one shown in Fig. |l|(a), with an array of cylindrical islands of 
arbitrary cross section, such as the ones shown in Fig. |l](b) and Fig. 0(a), provided that we 
keep the same junction capacitance C = eS/And and array period a. 

To bring the interaction problem to an analytically calculable form, we can investigate the 
continuum limit in which the discrete periodic structure is replaced by a continuous dielectric 
medium (e.g., as shown in Fig. [l](c) ) . This approximation is valid when the characteristic 
length of electron-electron interaction within the array is much larger than the array period, 
i.e. m ^> 1. The effective dielectric constant e e f of the medium can be found from the 
requirement that the relation between the average electric field (E) = (1/a) J a E dz and 
average displacement (D) is the same in the dielectric model as in the array model: 

a AnaC 

For the particular case of the dielectric model with circular cross section of radius 



R = yS/n (Fig. 0(c)), we can find the electron-electron interaction energy L^(r) from 
the electrostatic potential <p(p, z) induced in the dielectric cylinder by a charge e located at 
z = 0: 

U d (r) =e<f>(0,z = r). (4) 

For the relatively large distances we are interested in, m ~ m ^> 1 (i.e. ma ^> R) the 
shape of the initial charge is not important, and it is natural to spread it uniformly over 
a thin disk [z — 0,p < R). The resulting boundary electrostatics problem can be readily 
solved by the standard Fourier integral expansion: 

. 0(0,0) + ^[f™I (kp)A(k) cos(kz)dk-±\z\] , p<R, 
0(0, z) = { ef (5) 

/ °° K (kp) B(k) cos(kz) dk , P>R, 
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where the last term in the upper line takes care of the boundary condition (Gauss's law) 
on the charged disk at z — 0. On the other hand, the boundary conditions at p = R 

(d(j)/dz\ + = d(f)/dz\- 1 d(j)/dp\ + = e e / d<p/dp\-) give 
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A(k) 



7re ef k 2 



T(im+ h(kR)K (kR) 



(6) 



K^kR) 

[For numerical calculation of the integrals in Eqs. (H), @, it is convenient to use the 
expansion \z\ = (2/7r) J °° cos(kz) k~ 2 dk, in order to cancel the divergence of A(k) at k — > 0.] 
Numerical integration yields the functions U d (r) as shown in Fig. [| At large distances, 
these functions approach the free space Coulomb interaction: 

r 



U d (r)^-, r-^oo, (7) 



while at small distances, these functions are linear: 



U d (r) U d (0) - eE(0)r, r^O, (8) 
E(0) 2e 
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The crossover between these two limits takes place at r ~ r a = R^fe^j. Quite unexpectedly, 
at r ~ r the potential U d (r) is higher than the asymptotic value (0), approaching it at 
r —>■ oo from above. 

If r Q 3> a, we can apply this result to the discrete array of a circular cross section. In 
terms of the number of islands, the crossover point is 



m = — = — JeTf = \ — • (9) 



If we define the stray capacitance per unit period of the array as 



Co = l (10) 



we can once again formally express m as in the tridiagonal model, m = \/-£r- 

The dielectric-model potential function U d can be approximated reasonably well by a 
simple phenomenological expression: 



tt / \ e 2 ( a .—Km. 1 L ,-Km,ll .... 
U a (m) = -\—exp( ) + - 1-exp ) . (11) 

with two dimensionless fitting parameters, a and k. The boundary condition (J8|), together 
with the definition for m Q given in Eq. (|9|), yields the condition a = 2/k — k/2, leaving only 
one free parameter, k. Figure || shows the dependence of k on m , when k is adjusted to 
provide the best fitting (by sight) of U a {m) to Ud(m) (see Fig. |3|). One can see that k is 
close to unity, and depends on m only logarithmically, so that 

*»(« + *) „2e» /q_ e 2 

Thus, f/ a (0) may be approximated by Z7t(0) (Eq. (§)) if m Q is within the practical interval 
~ 3 — 10 (see below). 

In order to check the validity of our results for structures with a different cross section, 
we have calculated the interaction energy U s (m) for a chain of square cross section islands 
(Fig. 0(a)) for several values of the a/d ratio using the geometric capacitance modeling 
program FASTCApi. This program takes, as input, a collection of "panels," finite elements 
representing the surfaces of a group of conductors. By calculating the amount of charge 
induced on each panel when one conductor is held at fixed potential, the rest at zero, 
FASTCAP calculates, one by one, each row of the capacitance matrix for the group of 
conductors. 

For FASTCAP to calculate the capacitances accurately, the discretization of the surface 
must be fine enough to represent the charge distribution. The main criterion, when panels 
(and computer memory) are limited, is to let the panelling reflect the increased charge 
density near the corner of a conductor. Thus, when modeling conductors with square cross 
section, each face was divided into 9 panels (Fig. |](b)), with smaller panels along the edges 
where charge density is higher. To test the accuracy of this simple panelling, we modeled 
a chain of 10 islands in two ways: one with 9 panels per face, the other with 100 per face. 
Differences in electron-electron interaction energies between the two models were less than 
1%, leading us to believe that the 9-panel model is sufficiently accurate for our purposes in 
this work. 
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Island chains modeled with FASTCAP were limited in length to around 120 islands, due 
to computer memory limitations. To avoid edge effects, the energy of interaction between 
two electrons was calculated with both electrons separated from the edges of the array by at 
least ~ 2m islands. Energies calculated this way did not depend significantly on the total 
number of islands in the chain. The results of the FASTCAP calculations are shown as points 
in Fig. H One can see that they are very close to those obtained by the dielectric model, 
and hence to the analytical approximation of Eq. (|TT1) , despite a substantial difference in 
the geometry of the two cross sections. 

To summarize, we have shown that electron-electron interactions in long one-dimensional 
arrays of small tunnel junctions are much better described by Eq.flTTD (with k ~ 1, a ~ 2, 
and U a (0) ~ e 2 /2^/CC ) than by the traditional expression ([[]). The effective radius m = 
JcJC Q of the interaction is determined by the ratio of the tunnel junction capacitance C to 
the effective stray capacitance C Q fllCf). For an array formed on a substrate with dielectric 
constant e s , Eq. ([H]) should be modified to 

C * {es + 1)a . 
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For a typical present-day single-electronic array with a ~ 0.1/^m, formed on a SiC>2 
substrate (see, e.g., Ref. |7|), the formula above gives C Q /a ~ 10~ 16 F//^m and C a ~ 1CT 17 -F. 
Thus, for the typical junction capacitance C ~ 10 _16 F, the soliton radius m is close to 3. 
This estimate shows that for single-electronic devices where a considerable soliton radius is 
important (for example, for suppression of the macroscopic quantum tunnellingSH), vertical 
structures with stacked junctions may be more advantageous. Vertical stacking can provide 
a very small array period a, and hence the small stray capacitance C a needed for m ^> 1. 

The authors are grateful to J. White and K. Nabors for providing the program FASTCAP, 
and to D. Averin, A. Korotkov and J. Lukens for numerous discussions. This work was 
supported by AFOSR Grant No. 91-0445. 
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FIGURES 

FIG. 1. One-dimensional array of tunnel junctions: (a) a typical experimental geometry, with 
tunnel junctions of area S and permittivity e; (b) cylindric island model; and (c) dielectric approx- 
imation to cylindric islands. 

FIG. 2. (a) Square cross section island array model, (b) Boundary elements ("panels") for 
square cross section island model used for in FASTCAP capacitance calculations. 

FIG. 3. Functions £7(m) describing the shape of the single electron soliton. Solid lines: Ud(m), 
dielectric approximation to circular cross section islands, with e e f = 30,100,300 (bottom to top). 
To compare numerical and analytical results, we have set 7ri? 2 = a 2 , so that m = J e e j /vr. Squares: 
U s (m), square cross section array for the case a = b (see Fig. ||). Dashed lines: U a (m), as given 
by Eq. (|i~l~l). Dotted line: free space potential, U(m) = e 2 /ma. Note that a vertical offset is used 
in plotting the curves. 

FIG. 4. The best-fit parameter k, the sum k + a, and £7(0), as functions of m . Squares: 
£7 S (0). Open circles: £7^(0). Diamonds: £7 a (0). Error bars reflect the range over which a "good" 
fit, by sight, is achieved. 
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